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Abstract 



^ I The Minimal Supersymmetric Standard Model can possess several CP-violating phases be- 
yond the conventional Cabibo-Kobayashi-Maskawa phase. We calculate the contribution 
of these phases to T-violating nuclear forces. These forces induce a Schiff moment in the 
i99j^g nucleus, which is strongly limited by experiments aimed at the detection of the electric 
dipole moment of the mercury atom. The result for dng is found to be very sensitive to the 
CP-violating phases of the MSSM and the calculation carries far fewer QCD uncertainties 
than the corresponding calculation of the neutron EDM. In certain regions of the MSSM pa- 
rameter space, the limit from the mercury EDM is stronger than previous constraints based 
on either the neutron or electron EDMs. We present combined constraints from the mercury 
and electron EDMs to limit both CP-violating phases of the MSSM. We also present limits 
in mSUGRA models with unified gaugino and scalar masses at the GUT scale. 
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1 Introduction 



Despite the impressive success of the Standard Model, few are convinced that it is the 
final theory of particle interactions. For example, the supersymmetric modification of the 
Standard Model yields a very promising framework in which we are able to understand the 
stability of the electroweak scale. The Minimal Supersymmetric Standard Model (MSSM) 
provides a plethora of new phenomenological predictions which range from new charged and 
colored particles actively searched for in accelerators, to cold dark matter candidates, to new 
CP-violating phenomena such as the electric dipole moments of the neutron and electron 
which are generated if the additional CP-violating phases in the MSSM are non-zero. In this 
work, we study in detail the predictions of the MSSM for the electric dipole moment of the 
mercury atom and derive the constraints on the MSSM phases from the experimental limits 
on dHg- 

The null experimental results for the electric dipole moments (EDMs) of the electron, 
neutron, heavy atoms and diatomic molecules H, ^ Q can in general place very strong 
constraints on the CP-violating sector of a new theory and probe energy scales which are 
inaccessible for direct observations at colliders [^]. In general, the relevant contribution to 
the dipole moments at scales of ~1 GeV can be parameterized in terms of effective operators 
of different dimensions suppressed by corresponding powers of a high scale M where these 
operators were generated: 

A// = E ^.ot\ (1) 

n>4 

Here O'f''^ are operators of dimension n, with its field content, Lorentz structures, etc., 
denoted by i. The fields relevant for the low-energy dynamics of interest are gluons, the three 
light quarks, the electron, and the electromagnetic field. This general form is independent 
of the particular construction of the new theory, and the details of a given model enter only 
through the values of the coefficients c„j/M"~^. 

In the MSSM, the number of operators which can generate an EDM is considerably 
smaller than in the generic case. In fact, all four-fermion operators are numerically insignifi- 
cant. They can be generated in the MSSM only with additional factor of order {mq/ MsusyY 
modulo possible nontrivial flavor structure of the soft-breaking sector. Here we assume the 
minimal scenario with flavor-blind breaking of supersymmetry and therefore we can safely 
drop all four-fermion CP-violating operators. Hence, the relevant part of the effective La- 
grangian at the scale of 1 GeV contains the theta term, the three-gluon Weinberg operator, 
the EDMs of quarks and electron and the color EDMs (CEDMs) of quarks, 

i=u,d,s i=u,d,s 

We will assume here that the PQ mechanism of 6'-relaxation [^] eliminates 6 ~ 0(1) and 
sets 6 to Oeff at the minimum of the axion potential. When both the CEDMs and Weinberg 
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operator are absent, the value of 6'e// is exactly zero. However, nonzero w and di induce a 
linear term in the axion potential, and the effective value of 6 is different from zero. This 
value leads to an additional contribution to the EDM of the neutron, usually ignored in the 
literature. 

The coefficients in front of the operators in Eq. (H) can be calculated for any given model 
of CP-violation and then evolved down to the low-energy scale, using standard renormaliza- 
tion group techniques. In the MSSM, in particular, one can compute effective Lagrangian 
(§) for any given point in the supersymmetric parameter space. Then, to get the final predic- 
tions for EDMs, one has to take various matrix elements for these operators over hadronic, 
nuclear and atomic states 0, H, 0- In most cases this is a source of major uncertainty, 
especially when hadronic physics is involved. The exception is the case of a paramagnetic 
atom, in which the EDM is generated by the electron EDM de, and where the effects of 
nuclear CP-odd moments induced by the rest of the operators in can be safely neglected. 
The EDM of ^os^l is extremely sensitive to d^, due to a very large relativistic enhancement 
factor c ~ —600, which relates the EDM of the atom with de, d^i = cdg. The experimen- 
tal bound on the EDM of the thallium atom 0, combined with good stability of atomic 
calculations (see P] and references therein), leads to the following limit on the EDM of the 
electron: 

4 < 4 ■ 10-^^e ■ cm. (3) 

Therefore, the calculation of de in the MSSM gives the most reliable limits on CP- 
violating phases. It is clear, however, that the electron EDM limit alone cannot exclude 
the possibility of large CP- violating phases. This is because d^, as any other coefficient in 
Eq. (^, is in general a function of several CP-violating phases, and mutual cancellations 
are possible. This is what happens, for example, in the MSSM with the minimal number 
of parameters in the soft-breaking sector (see recent works |TI], |TT|). In the MSSM, it is 



well known that there are two independent CP- violating phases, 9^ and 9a, associated with 
the supersymmetric Higgs mass parameter /i and the soft supersymmetry breaking trilinear 
parameter A. The calculation of the relevant one loop diagrams determines de as a function 
of these two phases. If the phases are small, de is simply a linear combination of 9a and 
Therefore even a constraint as strong as that given in (^ leaves a band on the ^yi-^^ plane, 
along which a cancellation occurs and the phases are not constrained. In general, a second 
constraint could be expected to lift this degeneracy and place a strong constraint on both 
phases. It has been common to use the limit on the neutron EDM as this second constraint. 
Although there are large uncertainties in the calculation of the neutron EDM, as we argue 
below, when the limit on the neutron EDM is used, cancellations in the electron EDM occur 
in many of the same regions as cancellations in the neutron EDM. Therefore, one is led to 
the conclusion that large phases are still possible. 

In what follows, we critically reexamine the reliability of the calculation of the EDM of 
the neutron in the MSSM. We demonstrate that this calculation is subject to very large 
hadronic uncertainties, which makes the extraction of the limits on CP-violating phases in 
MSSM tenuous. Instead, we propose that useful limits may be obtained from the limits 
on the EDM of the mercury atom. This arises from the T-odd nucleon-nucleon interaction 
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in the MSSM, induced mainly due to the CEDMs of hght quarks. This interaction gives 
rise to an EDM of the mercury atom by inducing the Schiff moment of mercury nucleus. 
We demonstrate that the degree of QCD uncertainties related to this calculation is in fact 
smaller than in the case of dn and that it is possible to calculate the T-odd nucleon-nucleon 
interaction as a function of the different MSSM phases. As an example, we proceed with 
the calculation of the EDM of the mercury atom in one specific point of the supersymmetric 
parameter space where all squark, gaugino masses, \fi\ and |y4| parameters are set equal. This 
"pilot" calculation demonstrates the sensitivity of dng to the CP-violating phases of MSSM. 
We find in this case that dng provides somewhat better limits on CP-violating phases than 
dg. We proceed further and combine mercury EDM and electron EDM constraints to exclude 
most of the parameter space on 6 a- 0^ plane in this toy example. Finally, we consider more 
realistic constraints over the supersymmetric plane when supersymmetry breaking scalar and 
gaugino masses are unified at the GUT scale. In this case, we find that the limits on CP- 
violating phases obtained from dng is no longer more restrictive than (ig, as the RG evolution 
of soft-breaking parameters makes squarks and gluino heavier than sleptons, charginos and 
neutralinos. The combined limits are still very powerful as the cancellation of different 
supersymmetric contributions typically occur in different regions of parameter space. 



2 The Neutron EDM in the MSSM 

Limits on the neutron EDM are commonly used to set constraints on new CP violating 
interactions. In particular, the upper limit to dn is often used to limit the size of the CP- 



violating phases in the MSSM [T^. The current experimental limit on the EDM of the 
neutron is 

(i„ < 1.1 ■ 10-2^e- cm, (4) 

Indeed, the EDM of the neutron receives contribution from all operators listed above in Eq. 
(0) except rfg. However, there is a complication in using the neutron EDM as compared 
to the electron EDM, due to QCD uncertainties which make the extraction of the limits 
on CP-violating phases in the fundamental Lagrangian problematic. We demonstrate two 
aspects of this problem below. 

The most straightforward contribution to the EDM of the neutron is due to the quark 
EDM operators. It is usually estimated using nonrelativistic SU(6) quark model. The result, 

4 1 

dn — 2*^'^ '^^ui (5) 



can be compared, in fact, with the model calculations |13| and lattice simulations of light 



quark tensor charges in the nucleon [Q. The matrix elements for the tensor charges of the 
nucleon are defined by 

{N\i,,a^,'iP,\N) = 5qNa^,N, (6) 
whereas the axial charges are defined by 

(iVlV'.TMTs^.liV) = AgiV ^,^,N (7) 
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In the Naive quark model, both Lorentz structures correspond to the spin of a nonrelativistic 
quark. In this case 5u = Au = —1/3, 6d = Ad = 4/3, 6s = As = 0, yielding eq. (|^). 
Note that isospin symmetry gives us (An)„ = {Ad)p, etc. However, as argued in [|15|, 
since it appears that the contribution to the nucleon spin from the strange quark (A^) 



is no n- vanishing [16|, the naive quark model may not be sufficient to describe the quark 
EDM contribution to the neutron EDM. While it is not the axial charges which need to 
be considered for the calculation of the neutron EDM, but rather the tensor charges, the 
departure of the axial charge values from their NQM values indicates that more realistic 
(non-NQM) values of the tensor charges {6q) must be used. According to calculations based 



on Lattice QCD [|14|, the tensorial charges of up and down quarks in the proton should 
be read as Su ~ 0.8 and 6d ~ —0.23. This means that the naive nonrelativistic formula 
predicts the EDM of the neutron due to the quark EDMs to be 1.5-1.7 times larger than 
the lattice result. Slightly different values of ~ 1.1 and 6d ~ —0.4 can be derived from 
the SU(3) chiral quark soliton model |13[. The tensor charge of the strange quark is found 
to be consistent with zero in both methods |13, 14|. This is due to the fact that the sa^^s 



operator is odd under charge conjugation which must result in the Zweig-type suppression 
of this matrix element over the neutron state. Even with the usual rris/md enhancement 
of this operator, it is unlikely to be important. This does not exclude other possible CP- 
violating operators involving the s-quark, CEDM or generic four-fermion operators, as their 
contributions to the EDM of the neutron can be significant 0, 17|. Departures from the 



predictions of the non-relativistic quark model were recently considered in |jT8| . 

Unfortunately, the quantitative evaluation of the remaining contributions to the neutron 
EDM is complicated due to of our lack of knowledge about strong interaction dynamics 
at 1 GeV and below. Typically, one resorts to Naive Dimensional Analysis (NDA) [|19| , 
formulated within the constituent quark framework. This method is, however, only an order 
of magnitude estimate to be used when other methods of calculation fail to produce an 
answer. When the problem of estimating d^ due to du,d is considered, there are several 
possible answers in the literature: 



In 



e0.7{du + dd) Ref.[9] 



dN-^iO{l)d^ + 0{l)dd) NDA,Ref.[gg (9) 

We have chosen a normalization where gg is included in the definition of the operator in (0) 
and correspondingly include an additional Qs in the estimate @. The first result is based 
on a combination of chiral perturbation theory and QCD sum rules. The latter estimate is 
derived with the use of NDA For a realistic choice of the strong coupling constant at the 
scale of 1 GeV, Qs ~ VO.5 ■ Att = 2.5, the overall numerical coefficient in eq. is about 3.6 
times smaller than in (P). Estimates based on NDA imply that for natural relations among 
coefficients di/e ~ di, the effects of color EDMs on the electric dipole moment of the neutron 
are negligible and the result can indeed be approximated by the linear combination of EDMS 
of quarks. 



^We note that the estimate in ||2l|] is suppressed by an additional factor of 5s/47r. 
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In fact, it is possible to show that the CEDMs can lead to a substantially larger contri- 
bution to the neutron EDM than some of the predictions based on NDA. The easiest way to 
see that CEDMs can be numerically important is to calculate the effective 6'-term induced by 
CEDMs in the presence of the PQ symmetry and then use the result for dN{6). This value, 
6eff{di) can be calculated within the current algebra approach, in a manner similar to the 
calculation of the vacuum topological susceptibility p^, The dynamically induced 

theta term can be expressed in the following compact form: 



^e// = -^ f^ + f^ + f^ • (10) 




Here, is the ratio of the quark-gluon condensate to the quark condensate. It is known to 
good accuracy from QCD sum rules E^ that. 



2 _ {0\9qiGcr)q\0) 



^2 ^ ^ 1^-. y-^i / ^ o.SGeV^ (11) 



The accuracy of the estimate (|10D is of order m^^/m^/, which is acceptable for our 
purposes. If no interference with other terms is expected, then the expression (^) must be 
less than the current limit on 6, extracted from the same neutron EDM data. Using the fact 
that in the simplest variant of the MSSM, dd/rrid = ds/rris, and assuming for a moment that 
this is the only contribution to the EDM of the neutron, one can obtain the following, quite 
stringent, level of sensitivity for the CEDM: 

dd < IQ-^^cm. (12) 

This fact alone suggests that CEDMs may contribute significantly to the EDM of the neu- 
tron, typically at the level of the prediction (9) and an order of magnitude above NDA 



predictions. Remarkably, the main uncertainty in the limit (|T2D comes not from the cal- 
culation of 6{CEDM), but rather from the principal difficulties in calculating d]^{9). In 
the standard approach the chiral loop diagram is used to estimate d^iO). This loop is 
logarithmically divergent in the exact chiral limit and therefore is distinguished from the rest 
of the contributions. For realistic values of the parameters, however, this logarithm is not 
large and other contributions can be equally important. This makes the whole calculation 
problematic even in predicting the sign of the 9 term contribution to d^q. 

Besides dN{0{C EDM)), one should also consider direct CEDM-induced contributions to 
the EDM of the neutron which can be computed within the same chiral loop approach 0. 
Combining different contributions, we can symbolically write the result for the EDM of the 
neutron in the following form: 



dN ^ O.Sdrf - 0.23d„ + e 



du (ci In — + C2) + dd (c3 In — + C4) + d^ f C5 In + cq) 



(13) 

The coefficients Ci, C3 and C5 were estimated in Ref. to be Ci ln(mp/mjr) = C3 ln(mp/m^) ~ 
0.7 and C5 ~ 0.1. The cutoff parameter A corresponds to scales where chiral perturbation 
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theory breaks down, that is, of order rrip. In the exact chiral hmit, m^r, rriK 0, and the 
logarithmic terms dominate. In practice, however, the logarithmic terms are numerically 
not distinguished from the coefficients C2, C4 and ce, which are a priori comparable with 
Ci, C3 and C5 and are not calculable in this approach. It is clear then that these terms can 
change both the value and the signs of different contributions to d^- Therefore, although 
in principle very important as an order of magnitude estimate, Eq. (|13D fails to provide d]s[ 
as a known function of individual (ij-contributions and, ultimately, of different CP- violating 
phases. 



As emphasized in Ref. the NDA estimate of dniO) essentially reproduces the cal- 

culation of Ref. [^. The source of the disagreement in the case of dn{CEDMs) can be 
traced to the problem of estimating the CP-odd 7r"^pn-vertex, proportional to the matrix 
element {p\ugs{G(j)d\n) . In Ref. Q this matrix element was estimated to be -1.5 GeV^ and 
is essentially proportional to the quark-gluon condensate parameter ~ O.SGeV^ ([n]). On 
the other hand, it can be shown that NDA suggests for this matrix element a value of order 
i.e. one order of magnitude smaller. This difference is related to the 
fact that the NDA assumes nonrelativistic quarks whose chromomagnetic interactions are 
suppressed, whereas QCD sum rules use more realistic descriptions of hadronic properties in 
terms of vacuum quark-gluonic condensates. 

To summarize this discussion, the extraction of reliable limits on the CP-violating phases 
in the MSSM from the EDM of the neutron is difficult and uncertain. Even the best estimates 
of dn based on the "chiral logarithm" approach 0, bear a large degree of uncertainty and 
cannot produce a precise prediction for dn as a function of the CP-violating phases. Useful 
limits are still available from the electron EDM; however, the magnitude of the phases 
is not terribly constrained on this basis alone, due to cancellations in the various MSSM 
contributions to de- Fortunately, the EDM of the neutron is not the only source of information 
about CP-violation in the strongly-interacting sector. Limits on T-violating nuclear forces 
are provided by experiments aimed at the detection of the EDM of paramagnetic atoms, 
among which the EDM of ^^^Hg atom is the most constraining. In what follows, we will 
discuss the constraints these limits provide, both alone and in conjunction with the electron 
EDM hmits. 



3 CP-violating nucleon-nucleon interaction in MSSM 

The limits on T-odd nuclear forces extracted from the atomic experiments are in general very 
important for particle physics [0]. In the case of diamagnetic atoms, the most impressive 
hmit is obtained for the EDM of ^^^Hg 0: 

dHg<9- IQ-^^e ■ cm. (14) 

The electric screening of the electric dipole moments of the atom's constituents is violated 
by the finite size of the nucleus and can be conveniently expressed by the Schiff moment 
S, which parametrizes the effective interaction between the electron and nucleus of spin I, 
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Veff = —eS{l'V)6{r) [Q. Atomic calculations derive the atomic EDM as a function of S and 
translate the experimental result (|l^ into the limit on the Schiff moment of the nucleus: 



dHg = S- 3.2 ■ 10-i^fm~2 

5 < 2.8 • 10~^°efml (15) 

The Schiff moment of the nucleus can be induced either due to the Schiff moment of the 
valence nucleons or due to the breaking of time invariance in the nucleon-nucleon interaction, 
the latter being enhanced by the collective effects in the nucleus. The calculation of the Schiff 
moment of the nucleus, originating from various NNN'i'j^N' interactions was done in the 



single particle approximation with square- well and Woods-Saxon potentials |^8[. The results 
show that the Schiff moment of mercury is primarily sensitive to the ppni'^^n interaction. If 
we parameterize the coefficient in front of this interaction as ^Gp/ V^, the nuclear calculation 
p8| provides us with the following value for S: 



5 = -1.8-10"^ee-fm^ (16) 

Combined with Eq. ([T3|), it gives the following constraint on 

^<1.9-10~l (17) 

Questions concerning the calculation of the strength of ppni'-f^n interaction induced by dif- 
ferent operators was considered in @, ||, |^. The effective theta term, the Weinberg 
three-gluon operator and the CEDMs of quarks can generate this interaction. Numerically, 
the contributions provided by the CEDMs of up and down quarks are the most important 
and we concentrate our analysis on them, trying to incorporate the effect of dg as well. 



n 



n 



71 



-Tlpp ^'IPP 

Figure 1: Pseudoscalar meson exchange diagrams, inducing ppni'j^n interaction. 



Following [^, 0, we approximate the T- violating nucleon-nucleon interaction by pseu- 
doscalar exchange, as shown in Fig 1. In the limit of exact chiral symmetry this exchange 
has the power-like singularity to be compared with the logarithmic singularity in the 
case of the EDM of the neutron. The CP violation resides in proton-meson vertex which 
can be calculated with QCD sum rules and current algebra techniques. The CP-conserving 
meson-neutron vertex is sufficiently well known from SU(3)-relations in baryon octet decay 
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amplitudes and from the axial charges of nucleons. If only du and dd are present, the pion 
exchange dominates t] exchange by a factor m^/m^ ~ 16. In the MSSM, though, the strange 
quark CEDM is enhanced relative to that of the down quark by a factor irts/md and t] meson 
exchange is not a priori negligible. In the chiral approach, CP-violating vertices of interest 
can be reduced to the following set of matrix elements: 

Qirpp = ^" ^'^ (^{p\ugsiGa)u - dgs{Ga)d\p)) + 
^" ^'^ (^{p\ugs{Ga)u + dgs{Ga)d\p) — m'^{p\uu + dd\p) 



9vPP = {{p\sgs{.G(7)s\p) - rm? {p\ss\p)') (18) 

Here is the ratio of quark-gluon condensate to quark condensate introduced earlier in Eqs. 
(p!0|) and (11). At this point our results are already slightly different from 0, ^, 2^. Namely, 



we have included additional contributions related to the fact that the octet combination of 
color EDM operators has the quantum numbers of the 7r° and rj fields which can therefore 
be produced from the vacuum, vr", for example, can be "rescattered" on the nucleon with 
an amplitude proportional to (m^ + mu){N\uu + dd\N). As a result, the diagram shown in 
Fig. 2 is responsible for a contribution directly proportional to which is effectively of the 
same order as the direct contribution considered in H, 01. 



X 



N 




N 




<^oo 



Figure 2: Additional contribution to the g-i^NN, proportional to the nucleon sigma term 



Further calculation relies on QCD sum rules p2[ and low-energy theorems in QCD. 
Matrix elements from qgs{Ga)q operators were evaluated in ||^, |^, 



{p\qgs{Gcr)q\p) ^ -m^ {p\qq\p) ■ 



(19) 



The matrix elements over the proton can be obtained from baryon mass splittings and pion- 
nucleon scattering data. Here we take the following values for the matrix elements of qq over 



the nucleon 29 



{p\uu\p) ~ 4.8; {p\dd\p) ~ 4.1; {p\ss\p) ~ 2.8 



(20) 
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These values of qq matrix elements correspond to the choice = 4.5 MeV, = 9.5 MeV 
and TTis = 175 MeV. The values of these matrix elements, together with the factorization 
formula (p!9D, suggest that T-odd nucleon-nucleon forces are primarily sensitive to du — dd 
and insensitive to du + dd, simply because the contribution to g-j^pp proportional to du + dd 
in ([T8|) is relatively suppressed by 

9^pp{du + dd) 2{p\uu - dd\p) 

= — — _ J ~ U.z (zi) 

9npp{du - dd) {p\uu + dd\p) 

In this sense, the contribution furnished by 6eff is numerically insignificant b ecause Qnpp 
generated by 9 is also proportional to {p\uu — dd\p). 

Thus, these simple considerations suggest that due to the numerical dominance of the 
triplet combination of color EDM operators, the final answer for ^ takes the following form: 

e = G-/-^^{dd - du - 0.0124), (22) 

We can see that the contribution of the strange quark CEDM is numerically suppressed, 
mainly due to the additional smallness of rjNN CP-conserving interaction as compared to 

Combining equations (|l^), ( |16|) and (^2]), we arrive at the following prediction for the 



EDM of the mercury atom: 

dffg = -{dd -du- 0.0124) X 3.2 ■ lO^^e, (23) 

where the the numerical coefficient 3.2 ■ 10~^ corresponds to the choice of light quark masses 
given above. Using the experimental limits (p!^, we deduce the very strong constraint on 
the following combinations of the CEDMs of quarks: 

\dd -du- 0.01241 < 3.0 ■ lO"^'' cm. (24) 



It is important to note that the quark EDM operators cannot induce a large value for 5*. 
They do not induce the ni'-ffnpp interaction, and their contribution to the Schiff moment of 
the nucleus is associated only with electric dipole moment of the external valence nucleon 
0]. Current limits on dng are only sensitive to quark EDMs larger than 10~^^e ■ cm and thus 
these operators can be safely neglected. Similarly, the potential contribution from the three 
gluon operator GGG to dfjg are small. We rely here on the QCD sum rule estimates [pO| , 
showing no significant contribution from this operator to the T-odd nucleon-nucleon forces 
and thus to the EDM of mercury. 

Finally, we would like to comment on the accuracy of the predictions (^3]) and (|2^) , 
distinguishing between the error in the overall coefficients and the errors in the relative 
coefficients of 4-proportional contributions. The uncertainties of the atomic calculations of 
dug^S) and nuclear calculations of S{S,) mostly affects the overall coefficients. Although the 
uncertainty in the overall coefficient can be significant , it is acceptable for our purpose as 
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it influences only the width of the allowed region in 6^ — ^^-plane. What is more important, 
however, is that the prediction of the relative coefficients in front of individual di in eqs. (^) 
and (^) can be done in a more reliable way and we estimate that the accuracy of keeping 
the triplet combination dd — du and neglecting dd + du, eq. (plf) , is at the level of 20%. 
In effect, it makes the constraints imposed by dng much more useful than those provided 
by dn- Another advantage of the approach for calculating dng and dn, developed in refs. 
P, ^ and applied here, is that it reduces the error from the poor knowledge of the light 
quark masses. Indeed, even in the case of the naive formula for the EDM of the neutron, 
dn — {4:dd — du)/3, the individual quark EDM contributions are proportional to niu^d which 
are known to 50%. In the present approach, the answer for ^ is ultimately proportional to a 
linear combination of mj(0|gg|0), which can be rewritten as f^m^. times the function which 
depends only on the ratio of light quark masses, known to much better accuracy than the 
masses themselves. 



4 The limits on the MSSM CP-violating phases 



In previous work |Tl]], limits on the neutron and electron dipole moments were used to 
constrain the two independent phases (of fi and A) in the MSSM assuming that all the 
terms in the Higgs potential and all gaugino masses are real and that all of the A-parameters 
are equal at the GUT scale and share a common phase. In absolute terms, the phases 
are not overly constrained, ^ 0.3, for 6a — 7r/2. The reason for the lax limits, are 
several cancellations in the various contributions to the EDMs. Furthermore, in some regions 
of parameter space, these cancellations occur simultaneously for the electron and neutron 
EDMs. 

As we argued above, there are several reasons to suspect that the limit due to the neutron 
EDM must be treated with caution. Instead, we have argued that the limit coming from the 
EDM of Hg is the result of a "cleaner" calculation and carries fewer QCD uncertainties. In 
what follows, we will explore in detail the limits on the two phases using the constraint based 
on the EDM of ^^^Hg (|^) derived above. We will compare these constraints on the phases 
to that obtained from the electron EDM. As we will see, the cancellations in the EDMs do 
not always occur at the same points in parameter space. To demonstrate the importance of 
the mercury EDM limit, we first consider a SUSY model with a single mass scale. We then 
present general results which assume gaugino and sfermion mass universality at the GUT 
scale. 



Following ||3T| , p2| , we analyze the limits on 9a and 9^ for different values of super- 
symmetric parameters. To demonstrate the sensitivity of the mercury EDM to a com- 
mon scale of the supersymmetric masses with arbitrary and uncorrelated phases, we choose 
mj ~ Ma- ~ ~ \ Ak\ at the electroweak scale and take tan /? = 2. In Figures 3a and 3b, we 
show the sensitivity to the EDM of the mercury atom for the cases of sin^A = 1, sin^^ = 
and sin^A = 0, sin^^ = 1. At this particular point of the supersymmetric parameter space 
all of the calculations are significantly simplified. When all soft-breaking parameters are 
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sufficiently heavy, close to the TeV scale, the chargino and gluino propagators can be simply 
expanded m. vi/ M oi V2/M and only the zeroth and ffist order terms in the expansions need 
be kept. If needed, for lower values of gaugino masses, the results can be generalized to 
include all effects of mixing in the gluino and chargino sectors. 

The calculation of the chromoelectric dipole moments of quarks in MSSM was done in 



the series of papers pO|, |35|. When the CEDMs of quarks are induced by 9a, (as in Fig 
3a), the result is dominated by gluino exchange, with very small corrections coming from 
Ai-exchange: 

7 mrflAlsin^A / ^gl gl \ 

J ^ _ mu\A\ sin Oa + ^\ 
^ 167r2M3 1^ 18 54/ 

Here rj denotes the renormalization group factor which reflects the QCD evolution of the 
color EDM from the weak scale to 1 GeV. When the color EDM operator is defined as in 
eq. (1^), its anomalous dimension is negative and small so that the overall renormalization of 
di is not important. The alternative definition of color EDM operator, frequently occurring 
in literature, is \d' qf^G^^j^a ^j.^'-^^q, where Qs is included in d' . Defined this way, this operator 
acquires a renormalization factor roughly proportional to g s{\-GeV) / gs{Mz) — 2 which is 
smaller than the value 3.3 quoted in [Q. This is because in Refs. ||20|, |3^ a very large 
coupling constant at low energies, — 2, is used. There is, however, an important numerical 
contribution to 77 which reflects the suppression of light quark masses at the high energy scale, 
mci{Mz) I niciilGeV) . We choose to use low energy values for m„ and m^, 4.5 and 9.5 MeV, 
and include the quark mass RG factor into rj. For the scale M of order M^, rj is numerically 
close to 0.35 and is mainly due to the suppression of the quark masses at the high energy 



scale. This suppression factor was omitted in Ref. [|20| where mu{Mz) = 8 MeV is used. 
Combining all numerical factors, we obtain the following value for the EDM of ^^^Hg: 

d^, = e. 1.5 . ^Q-^^ irn.-rn.-0M2m.)\A\sineA ^ ^ . j^mV^|^ ^ . ^^6) 

where we simply take \A\ = M. We see from Fig. 3a that the mercury EDM places a 
constraint on M, M > 1.5 TeV. 

In the other case, with sin^^ = 0, sin^^^ = 1, we have to include A2-higgsino and Ai- 
higgsino exchanges as well, so that the result for the CEDMs, (shown in Fig. 3b), is as 
follows: 

~ md|/i| tan/5sin6'„ , w., ^1 , , 

^'^^^ Lm^ H t^ + V + T^I (27) 




~ _ m„|/i|cot/3sing^ /5^| g^ 7gf\ 
" 167r2M3 1, 18 ^ 8 ^ 216 j ■ 

As a result, the contribution of the up quark relative to that of the down quark is suppressed 
by m„/ (m^ tan^ (3) . Numerically, the gluino exchange diagram dominates again with less than 
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Figure 3: The sensitivity of the EDMs of mercury and electron to the scale of the soft- 



breaking parameters with a) maximal phase of A [6 a = 7r/2, 9^ 
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a 10% contribution coming from A2-higgsino exchange. In this case, the hmit is somewhat 
stronger giving, M > 3 TeV. 

As one can see, the EDM of mercury is sensitive to the scale of supersymmetric masses 
as high as 1.5-3 TeV. This can be compared with the sensitivity of the EDM of the electron, 
which we calculate at the same point of the supersymmetric parameter space, taking the 
slepton masses equal to the squark masses: 



de 



dp 



m,.\A\ sin ^.4 gf 
167r2M3 12 
ruelnl tan /9 sin 9^ 



(28) 



24 24 



The limits based on the electron EDM, for the two cases considered, are weaker as can be 
seen from Fig. 3a and 3b where the limit on M is 0.4 and 1.7 TeV. 

There is also the possibility of destructive interference between two contributions induced 
by the CP-violating phases. Again, we choose the supersymmetric parameters to be equal 
and fix them to be in the range from 250 - 750 GeV. Fig. 4a-4c show the combined exclusion 
plots. The two bands correspond to the parts of the parameter space where the mercury 
or electron (Tl) constraints are hfted by the cancellation of different supersymmetric con- 
tributions. The allowed area lies on the intersection of these two bands. We observe that 
the band corresponding to the mercury EDM constraint has a different slope than that of 
the electron EDM, mainly because dng is by far more sensitive to 9a- We observe that both 
phases are sufficiently constrained for the low values of M. 
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Fig. 4a: M =250 GeV 



Fig. 4b: M =500 GeV 
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-0.1 -0 . 05 0.05 0.1 

Fig. 4c: M =750 GeV 

Figure 4: Combined, dug and c/e, constraints on the supersymmetric phases Oa/t^ and ^^/vr 
for different values of M. Allowed area is on the intersection of two bands. 

5 EDMs in mSUGRA and Cosmological Constraints 

We now consider the constraints on CP violating phases in mSUGRA-like models, i.e. models 
with unified gaugino and sfermion masses. We recall that to one loop, the phase of /i does not 
evolve with scale, but the phases of A^t and A^. must be run separately from the unification 



scale to low energies. We follow the analysis of [|TT], but with two changes. First, we 
replace constraints from the neutron electric dipole moment with limits from the EDM of 
Hg, discussed above. Second, we include recent results on the effect of coannihilations of 
neutralinos with staus on the neutralino relic density [Q. The latter has the effect of 
weakening the cosmological upper bound on the gaugino masses. This is demonstrated in 
Fig. 1^, where the light shading indicates the region of the {mo, "021/2} plane which yields a 
neutralino relic abundance in the cosmologically preferred range 0.1 < ^^h"^ < 0.3. The 
upper limit of the light shaded region crosses below the line = at mi/2 ~ 1400 GeV; 
for greater mi/2, either the relic density violates the upper bound ^^h"^ < 0.3 (which follows 
from a lower limit of 12 Gyr on the age of the universe) or the lightest supersymmetric particle 
is a stau, leading to an unacceptable abundance of charged dark matter. Here we've taken 
tan/3 = 2, but the light shaded region is quite insensitive to tan/3 for the values of tan/3 



we consider, as well as insensitive to the phase of For comparison, the dashed lines 

demarcate the inferred cosmologically preferred region if one ignores the effects of neutralino- 
slepton coannihilation. Whereas in [11|, the constraint ^^^^ < 0.3 yielded an upper bound 
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of 450 GeV on mi/2, we now have to consider larger values of mi/2- However, we will see that 
this does not effect the upper bound on 6^. 



200 i ; ' ! ; :' 



o 




100 200 300 400 500 600 700 800 

Figure 5: The light-shaded area is the cosmologically preferred region with 0.1 < VL~ h"^ < 0.3. 
The dashed line shows the location of the cosmologically preferred region if one ignores the 
light sleptons. In the dark shaded region the LSP is the f^, leading to an unacceptable 
abundance of charged dark matter. Also shown as a dotted line is the contour m^± = 95 GeV. 



In contrast to the results of the previous section, we find that in mSUGRA-like models, 
constraints from the electron EDM are typically more restrictive than those from the EDM 
of Hg. This difference arises because in models with gaugino masses unified at the GUT 
scale, the gluino tends to be considerably heavier than the neutralino and charginos, and 
this suppresses the contribution to dng from the quark chromoelectric dipole moments due to 
gluino exchange. We recall that cancellations between the chargino and neutralino exchange 



contributions to the electron EDM allow for large values of ||32| , |35| , |TT[. A similar effect 
also applies in the case of the Hg EDM, where cancellations can occur between the gluino 
exchange and neutralino and chargino exchange contributions to the quark chromoelectric 
dipole moments. The power of combining the electron and Hg limits lies in the fact that for 
fixed 9^ and 9 a-, the cancellations in the electron and Hg dipole moments occur for different, 
and often non-overlapping, ranges in my 2 ■ Thus the combined limits are stronger than either 
limit alone. 



Following ||32| , |TT[, we compute the electron and Hg EDMs in mSUGRA as a function of 
9f^,9A and mi/2 for fixed AQ,mo and tan/3. In Fig. ^-c we display the minimum value of 
mi/2 required to bring both the electron and Hg EDMs below their experimental limits, for 
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tan (3 = 2 and ttlq = 130 GeV. We exclude points which violate the current LEP2 chargino 
and slepton mass bounds fS^. The EDMs are computed on a 40x40 grid in {^^,6*^}, and 
features smaller than the grid size are not significant. Although the dependence of the EDMs 
on mi/2 is not monotonic, there is still a minimum value of mi/2, due to cancellations, which 
is permitted. In the zone labeled "1", m"^^ < 200 GeV, while the zones labeled "11", "111", 
"IV" and "V" correspond to 200 GeV < m"^^ < 300 GeV, 300 GeV < m"^^ < 450 GeV, 
450 GeV < m'lji^' < 600 GeV and m'^i^ > 600 GeV, respectively. Comparing with Fig. ||, 
we see that values of mi/2 larger than about 600 GeV are cosmologically excluded for this 
value mo. Therefore, region V corresponds to an excluded region in the phase plane. Of 
course, for this value of tan/?, the current Higgs mass bound requires enormous sfermion 
masses ^ 1 TeV, which are cosmologically prohibited. We've chosen to plot our results 



for tan/? = 2 in order to compare to our previous results p2|, [TTJ. Qualitatively similar 



conclusions apply for larger tan /?, which we summarize at the end of this section. 



Figure 2 of Ref. [jTT| displays the corresponding contours to our Fig. ya-c, but imposing 
only the constraint from the electron EDM|. Note that in |ri| we do not include a contour 
corresponding to 450 GeV < m™™ < 600 GeV, as this region would be cosmologically ex- 
cluded in the absence of coannihilations of neutralinos with sleptons, whose effects were not 



included in [|rT|. The effect of including the Hg EDM bounds is particularly significant at 
large Aq, where the cancellations are enhanced and the bounds on 6^ are weakest. Here the 
widths of the allowed region in my 2 at fixed 6a and 6'^ are narrowest, leaving less oppor- 
tunity for overlap between the ranges allowed by the electron and Hg EDMs, respectively. 
Indeed, for Aq = 1.5 TeV, the upper bound on 6*^ is reduced from ~ 0.37r, in the case of 
the electron EDM alone, to ~ O.lSvr, combining the two constraints, and, further, the width 
of the region in is considerably narrowed. The reduction in the bound on 6^ is minimal 
for small Aq, where the bounds on 6'^ are strongest. However, notice that the m™™ at the 
largest allowed values of 6fj_ is shifted from less than 200 GeV in the case of the electron EDM 
alone to between 200 and 300 GeV for the combined bound, in the case Aq = 300 GeV. For 
Ao = l TeV and 1.5 TeV, m^^J^'^ lies above 300 GeV at the largest O^,. 

We note that the larger values of mi/2 which neutralino-slepton coannihilation permit 
do not increase the maximum 6'^, for this value of mo. This is because, as we see from 
Fig. H, the region of mutual cancellations happens to lie at lower mi/2, between 300 and 
400 GeV. The widths of the allowed regions in mi/2 are typically between 50 and 80 GeV for 
the lightest shaded zones in Fig. and ^c and greater than 80 GeV almost everywhere in 
Fig. Larger mi/2 does, however, widen the allowed swath in by the region labeled 
"IV" in Fig. 1^. It helps in particular at small where the electron EDM can be beaten 
down sufficiently by taking heavy gaugino masses, without resorting to cancellations between 
different contributions. At large mi/2, the Hg EDM typically provides little constraint on 6^ 
due to the heaviness of the gluinos. As mo is increased, the regions of cancellation shift, and 
the maximum value of 6'^ slowly decreases. 

These effects are enhanced at larger mo and mi/2, in the cosmologically allowed "trunk" 



^In ||ri| we take toq = 100 GeV, rather than 130 GeV; however, taking mo ~ 130 GeV makes only a small 
change in the displayed contours and a slight reduction in the upper bound on 0^. 
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Figure 6: Contours of m™™, the minimum mi/2 required to bring both the electron and 
Hg EDMs below their respective experimental bounds, for tan/5 = 2, mo = 130 GeV, and 
a)ylo = 300 GeV, h)Ao = 1000 GeV and c)Ao = 1500 GeV. The central light zone la- 
beled "I" has m"^^ < 200 GeV, while the zones labeled "II", "III", and "IV" correspond 
to 200 GeV < m"^^ < 300 GeV, 300 GeV < m"^^ < 450 GeV, 450 GeV < m^^° < 600 GeV 
and m™™ > 600 GeV, respectively. Zone V is therefore cosmologically excluded. Panel d) 
shows the allowed region in the "trunk" for Aq = 300 GeV, mo = 200 GeV. The light shaded 
regions "A" and "B" are permitted, whereas the dark shaded region is cosmologically ex- 
cluded (see the text). 
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Figure 7: The maximum values of 6^ allowed by cosmology and both the electron and Hg 
EDMs, as a function of tan /3, for mo = 100 GeV (thick hues) and mo = 200 GeV (thin lines) 
and for Aq = 300, 1000 and 1500 GeV. 



which hes on top of the f/j LSP region (see Fig. The allowed region narrows as tuq 
increases, and the ^^h'^ = 0.3 contour crosses the line m^^ = and gives an upper bound 
on mi/2 and mo at mi/2 ~ 1400 GeV, rriQ ~ 300 GeV. The trunk region yields much larger 
sparticle masses than are cosmologically permitted in the absence of coannihilations, and this 
can suppress the contributions to the electric dipole moments sufficiently so that significant 
cancellations between the various contributions are not necessary. For low Aq, where the 
bounds on 6^ are tightest, the bounds on 6^ are somewhat relaxed in the trunk area. In 
Fig. ||d, we display the allowed region in the {9^, 9a} plane for Aq = 300 GeV, mo = 200 GeV. 
For this value of mo, my2 is cosmologically restricted to lie between 850 GeV and 950 GeV. In 
the light region labeled "A" , the EDMs are below the experimental limits for all 850 GeV < 
^1/2 < 950 GeV, while in the regions labeled "B" , only part of this range of my2 satisfies 
the EDM constraints. The dark regions at large \9^\ require my2 > 950 GeV to satisfy the 
EDM bounds, and so these regions are cosmologically excluded, as they yield a stau LSP. 
The upper bound on is relaxed to ~ 0.055. Taking mi/2 and mo at their maximal values 
allows 9fj_/n up to about 0.1. 

For large Aq, where 9^ can take its maximal values, the bound on does not weaken in 
the trunk region. As above, this is due to the fact that at larger 9^, cancellations are still 
required to bring the EDMs below their experimental limits, and the regions of cancellations 
occur at lower mi/2. Even taking mi/2 and mo at their largest cosmologically permitted 
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values does not allow for 9^ larger than the bounds in Fig. ^3,0. Further, since the regions 
of low my 2 are cosmologically forbidden at large mg, the bounds on 6^ at large Aq actually 
decrease for large rriQ. Thus the presence of the coannihilation trunk region does not increase 
the overall combined cosmology/EDM bound on the phase 9^. 

Lastly, we plot in Fig. ^ the maximum value of 9^ allowed by the electron and Hg electric 
dipole moments and the upper limit on as a function of tan/?. The thick lines are for 

mo = 100 GeV, while the thin lines are for mo = 200 GeV and show the effect on the bounds 
described above as one moves into the trunk region. 

6 Conclusions 

We have shown that the calculation of the EDM of the neutron as the function of different 
MSSM phases is problematic due to large uncertainties related to the contributions of the 
color EDMs. This is in contrast to the electric dipole moment of the mercury atom, induced 
by the T-odd nucleon-nucleon interaction. In the chiral limit the coefficient ^, characterizing 
the strength of T-odd forces has a power-like singularity ~ m'"^ ^ whereas logm^ in 

the same chiral approach. It is apparent that the vr'' and 77 exchange diagrams dominate 
both parametrically and numerically and therefore yield a very good approximation to the 
magnitude of the T-odd interaction. The final result is proportional to {da — d^ — 0.012ds) x 
3.2 ■ 10~^e and can be further developed in terms of CP-violating phases of the MSSM. 

There are two serious problems with the calculation of the T-odd nuclear forces due to 
the effective interaction (|^) with the coefficients provided by the MSSM. First is the status of 
the factorization in Eq. (^9|), related with the low-energy theorem in O"*" channel. Following 
Refs. 0, H, |], we have have taken {p\qgs{Ga)q\p) ~ 1.3GeV^(p|gg|p). We note that a 
designated sum rule calculation of this quantity and/or its simulation on lattice is highly 
desirable for it is the main source of uncertainties in the calculation of T-odd nuclear forces. 
The second potentially troublesome point is the effective negative sign between the da and 
ds contributions. Although the numerical suppression in front of dg is quite strong and dd 
dominates, destructive interference is still possible in both cases, 9a 7^ and 9^ 7^ 0. 

In this paper, we have considered first a very specific part of the supersymmetric pa- 
rameter space, when all squark, slepton and gaugino masses were chosen to be equal. The 
theoretical prediction for dng exhibits remarkable sensitivity to the scale of soft-breaking 
mass parameters as high as 1.5-3 TeV. When the scale is fixed below 1 TeV, dng limits both 
phases. The constraints on the CP-violating supersymmetric phases, obtained in this way 
are the strongest constraints so far. 

We have also considered the combined constraints from the Hg and electron EDMs in the 
mSUGRA, when all supersymmetry breaking gaugino masses, soft scalar masses, and soft 
trilinear terms are separately unified at the GUT scale. In this case, the sensitivity to the 
Hg EDM is weakened due to the relative size of the gluino mass. Nevertheless, the results 
are as strong or stronger (particularly when |y4| is large and the limits are weakest) than 
the combined results from de and d^- The improvement in the limit is due to the fact that 
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cancellations among the contributions to the EDMs occur at slightly different regions of the 
SUSY parameter space. 
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